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INTRODUCTION 



In any first order phase transition a quantity of fundamental importance is the rate 
at which bubbles of the equilibrium ( "true vacuum" ) phase nucleate within the metastable 
("false vacuum") phase. The bubble nucleation rate per unit volume, T, can be calculated 
^ ■ by a method, due to Coleman* which is based on finding a "bounce" solution to the 
classical Euclidean field equations. Thus, for a theory with a single scalar field, one seeks 
a solution of 

M M * = (1) 

for which the is near its true vacuum value near the origin, but approaches its false 
vacuum value as any of the tend to ±oo. 

A difficulty arises when one deals with theories where the symmetry breaking is a 
result of radiative corrections? In such cases the vacuum structure is not determined by 
V((f>), but instead can be found only by examining the effective potential, V e ff((f)). The 
bounce equation (1) is clearly inappropriate — in fact, if V (0) has only a single minimum 
there will not even be a bounce solution. An obvious solution is to replace V((f)) by V e s ((f)) 
in the bounce equation? 



Based on a talk given at the NATO Advanced Research Workshop Electroweak Physics 
and the Early Universe, Sintra, Portugal, March 25, 1994 
This work was supported in part by the US Department of Energy 



1 



Although plausible, and clearly a step in the right direction, this procedure raises 
some questions. The one-loop radiative corrections generate an effective action which 
contains not only V e g, but also terms involving derivatives (of all orders) of the fields. Can 
these terms be neglected when dealing with configurations, such as the bounce solution, 
which are not constant in (Euclidean) space-time? Even if this can be done in a first 
approximation, what are the nature and magnitude of the corrections which these terms 
generate? There are also questions relating to V e & itself. First, the effective potential 
obtained by perturbative calculations differs considerably from that defined by a Legendre 
transform (the latter must be convex, while the former is not). The latter clearly does not 
lead to an appropriate bounce, but how precisely does the formalism pick out the former? 
Further, the perturbative effective potential is known to be complex for certain values of 
the fields. How is the imaginary part of V e s to be handled? 

In this talk I will describe a systematic calculational scheme 4 which gives an answer 
to these questions. The general idea is to use the path integral approach of Callan and 
Colemanf but to integrate out certain fields at the outset. This leads to a modified 
effective action that gives a correct description of the vacuum structure of the theory 
and has a bounce solution that can provide the basis for a tunneling calculation. To 
leading approximation, one obtains the same result for V as would have been obtained by 
replacing V by V e g in the standard procedure. The next-to-leading terms, which can be 
expressed in terms of the effective action, give calculable and significant corrections. The 
corrections beyond these, although well-defined, do not have a simple expression in terms of 
the effective action. In particular, the potentially complex terms in the effective potential 
do not appear directly, but only as part of more complicated functional determinants that 
can easily be shown to be real. 

Throughout this talk I will confine myself to the case of bubble nucleation by quantum 
mechanical tunneling at zero temperature. Similar issues arise in connection with finite 
temperature bubble nucleation, even in theories where radiative corrections have little 
effect on the zero-temperature vacuum structure. It should be possible to deal with these 
by appropriate extensions of the methods described here. 

A SIMPLE MODEL 

Although radiative symmetry breaking is of greatest interest in the context of gauge 
theories, I will begin by using a somewhat simpler model to illustrate the method. This 
avoids some technical complications associated with gauge theories, which I discuss later. 
The model has two scalar fields, <j> and A, and is governed by the Lagrangian 

C = l(d,<P) 2 + ±(d,A) 2 -V(<p,A) (2) 

with 

V = + ±<f + \m*A> + LA* + \g^A\ (3) 

Both fj 2 and m 2 are positive, so the tree-level potential has only a single symmetric min- 
imum. We would like A-loop corrections to induce symmetry breaking with ((f)) = a and 
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Figure 1. The two-loop graphs which contribute to the effective potential at order g . Solid and dashed 
lines represent (j) and A propagators, respectively. 

(A) = 0. This is done by choosing parameters so that n = 0(g 2 a), m = O(ga), f = 0(g 2 ), 
andA = 0(# 4 ). 

The one-particle-irreducible Green's functions of this theory are generated by the 
effective action S e s, which can be expanded in the familiar derivative expansion 



SeffOM) = / d 4 X 



V e g((f), A) + l -Z^A){d^f + \z A {^A){d^A) 2 + 



(4) 



where the dots represent terms containing four or more derivatives. The leading con- 
tribution to the effective potential is of order g 4 . It is given by the sum of the tree- level 
potential, all graphs with a single A-loop, and the corresponding counterterms and is equal 
to 



I 



(V; ff ) 3 4(^) = --(m^ + ^ 2 ) 2 



64tt 2 



In 



m 2 + g 2 4> 2 
m 2 + g 2 a 2 



(5) 



where y\ is a constant of order g 2 m 2 /a 2 . The order g 6 contributions to V e ff arise from the 
two graphs shown in Fig. 1, together with several counterterm graphs. In these graphs 
the propagators are not simply those of the tree-level theory, but take into account the 
interaction with a background 4> field. Thus, the A propagator is given by 

G A (k 2 ) = i(k 2 -m 2 -g 2 ( p 2 )- 1 . (6) 

The 4> propagator is a bit more complicated. Because the tree-level and one A-loop con- 
tributions are of the same magnitude, the effects of both must be included, giving 

G4k 2 )=t[k 2 -v 2 -(V eS )^}- 1 . (7) 

(The inclusion of one-loop effects in this propagator means that when calculating higher 
order corrections one must be careful to avoid double-counting of graphs; this will not be 
relevant for the calculations discussed here.) At order g 8 there are both three- loop graphs 
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containing A propagators and graphs with a single (ft loop. The latter graphs bring in 
terms proportional to ln(l^,ff)" 4 that become complex for certain values of (ft. 

This complex effective potential deserves some comment. The effective potential is 
often defined formally through a Legendre transform of the generating functional of con- 
nected Green's functions. The properties of Legendre transforms then imply that it is 
everywhere convex, a property that the perturbatively calculated effective potential does 
not enjoy in theories with multiple stable or metastable vacua. To understand this, re- 
call that the effective potential V^ g (4>) obtained by Legendre transform is equal to the 
minimum value of the energy density among all states such that (^>\<ft>(x)\$>) = (ft. For 
values of (ft that lie between two vacua, the energy is minimized by states |^) which are 
superpositions of the two vacua; these lead to a flat effective potential. The perturbative 
effective potential is most easily understood 6 by decomposing </>(x, t) into a spatially uni- 
form mode 4>o(t) and a part 0(x, t) whose spatial integral vanishes. In the infinite volume 
limit the former mode is essentially classical, and one can discuss states in which the wave 
functional is of the form S((fto — 0)\I/[/i/>(x]. The real part of the perturbative effective po- 
tential is the minimum expectation value of J;he energy density among states of this form 
subject to the additional condition that ^[^(x] be concentrated near ift = 0. For values 
of (ft between two vacua these latter states are not eigenfunctions of the Hamiltonian and 
eventually decay to states whose wavefunctionals are concentrated on configurations with 
large fluctuations in ift(x.). Being unstable, they have complex energies with the imaginary 
part of the energy related to their decay rate. In the region where the tree-level potential 
is concave the decay is essentially classical and is reflected in a perturbative imaginary 
part to the effective potential. Beyond this region decay proceeds by quantum tunneling 
and leads a nonperturbative imaginary part. 



THE CALLAN-COLEMAN FORMALISM 

In theories where the vacuum structure can be read off from the tree-level Lagrangian, 
r can be calculated using a formula that was derived by Callan and Coleman 5 using path 
integral techniques. The starting point is the quantity 



G(T) = ^(x)=^ fv |e- HT |^(x) 
= J [d(ft}e-^ +s ^ 



(8) 



where the Euclidean action 



r r T/2 

S((ft) = / d s x / dx A 

J J -T/2 



l(d,(ft) 2 + V((ft) 



(9) 



and S ct contains the counterterms needed to make the theory finite. 

The path integral is over all configurations such that (ft takes its false vacuum value 
0f v at X4 = ±T/2 and at spatial infinity. In the limit T — > 00, Eq. (8) is dominated by the 
lowest energy state with a non- vanishing contribution (i.e., the false vacuum) and is of the 
form 

G(T) « e~ £Tn (10) 
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where O is the volume of space and 8 may be interpreted as the energy density of the 
false vacuum state. Because this is an unstable state, £ is complex with its imaginary part 
giving the decay rate, which in this case is simply the bubble nucleation rate. Dividing by 
O gives the nucleation rate per unit volume, 

r = -2Im£. (11) 

The path integral may be approximated as the sum of the contributions about all 
of the stationary (or quasi-stationary) points of the Euclidean action S '(0): the pure false 
vacuum, the bounce solution 4>b with all possible locations in Euclidean space-time, and all 
multibounce configurations. In each case the contribution to the path integral is obtained 
by expanding the field about the classical solution (p(x): 

4>(x) = 4>{x) +r](x) (12) 

and then integrating over rj. To leading approximation one keeps only the terms in the 
action which are quadratic in rj. Expanding these in terms of the normal modes of S"(4>) = 
— □ + V"{4>) gives a product of Gaussian integrals. The evaluation of these integrals is 
completely straightforward in the case of the false vacuum, but about the bounce solution 
is complicated by the fact that S"(4>b) has four zero and one negative eigenvalues. The 
zero modes are treated by introducing collective coordinates; integrating over these gives a 
factor of OT. The negative mode is handled by deforming the contour of integration; aside 
from a factor of 1/2, this gives a contribution whose imaginary part is just that which 
would have been obtained from a naive application of the Gaussian integration formula. 
Finally, the contributions from the multibounce configurations are simply related to that 
from the single bounce. Summing over all of these stationary points gives 

G = Go + nTG G b + \{QT) 2 GoGl + ... = G e nTGb . (13) 

This leads to 

T = 2\mG h = e - B KJe- [SctiM - Sct{M] (1 + • • •) (14) 

where 

B = 5(0 b ) - S(0 fv ) (15) 

and 

det'[-n+V"(<P h )} 

det[-n + y"(0 fv )] 

Here det' indicates that the translational zero- frequency modes are to be omitted when 
evaluating the determinant. The determinants are divergent, but these divergences are 
cancelled by the terms containing S c t- Finally, J is a Jacobean factor associated with the 
introduction of the collective coordinates; it can be shown to equal B 2 /An 2 . 

Scaling arguments can be used to estimate the magnitude of the various terms in the 
expression for V. Let us assume that we can identify a small coupling A and a dimensionful 
quantity a such that the potential can be written in the form 

V(<P) = Xa 4 U(i;) (17) 



(16) 
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where U involves no small couplings and the dimensionless field ip = <p/a. The minima of 
the potential must then be located either at = or at values of <f> of order a. 

By defining a dimensionless variable s = \f\ax, we may write the field equations as 

dU , s 

a* = w (is) 

From the assumptions made above, this equation involves no small parameters and so has 
a bounce solution in which if; is of order unity and differs from the false vacuum within a 
region of a spatial extent (measured in terms of s) which is also of order unity. In terms of 
the original variables, the bounce has (f> of order a and extends over a range of x of order 
l/(V\a). 

With the same change of variables, the action becomes 



dip 

a^ 7 



+ U(if>) 



(19) 



Since the integrand contains no small parameters, while the volume of the bounce restricts 
the integration to a region of order unity, the bounce action B is of order A -1 . 
Similarly, the determinant factor K becomes 



K = - 



det'[(A<7 2 )(-Q, + *7"(V>b))] 



det[(Aa 2 )(-D s + £/"(Vfv))] 
det'[-n s + U"(ip h )} 



-1/2 



1 \ 2 ^4 

2 Xa 



det[-n s + *7"frM] 



1/2 



(20) 



where the explicit factor of A 2 cr 4 on the second line arises because the det' factor involves 
four fewer modes than the det factor. With this factor extracted, the ratio of determinants 
is formally of order unity, although divergent. Finally, the Jacobean factor J is proportional 
to B 2 ~ A -2 . Putting all of these factors together, we see that the nucleation rate is of 
the form 

T = aa 4 e~ b / X (21) 

with a and b both of order unity. In practice b can be calculated rather accurately, since 
it is not difficult to numerically solve the differential equation for the bounce solution. On 
the other hand, an accurate determination of a, which involves a functional determinant, 
is exceedingly difficult. 



A MODIFIED FORMALISM 

If we were to apply this formalism to our scalar field model, the first step would be to 
find a bounce solution to the Euclidean field equations implied by the Lagrangian (2). The 
problem is that, because the tree level potential has only the symmetric minimum, there 
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is no bounce solution. We can circumvent this problem by integrating out the A field from 
the start, and writing 



G(T) = J [d<P][dA]e- s ^' A) = J [d(j)]e 



W(4>) 



(22) 



W may be thought of as a kind of effective action, although it is not the same as the more 
usual S'eff of Eq. (4) . Graphically, the two actions differ in that W receives contributions 
only from graphs with only internal A-lines and external 0-lines. Also, W is divergent 
while S e ff, since it generates renormalized Green's functions, must be finite. Performing 
the path integral over A fields, one obtains 



W(<j>) = J d 4 x 



+ ^ V + + \{x\ In [-□ + M' 



\x) 



(23) 



+ counterterms + 



where M 2 ((p) = m 2 +g 2 (p 2 and the dots represent the two- loop and higher order corrections. 

Although W ', in contrast with the tree-level action, does reflect the true vacuum 
structure of the theory, we cannot simply proceed by solving its field equations to obtain 
a bounce. There are two difficulties here. First, we only have a perturbative expansion 
for W. Second, and more importantly, W is nonlocal and will therefore lead to quite 
complicated field equations. 

Although nonlocal, W can be approximated by a local functional if <fi is sufficiently 
slowly varying. (For the terms shown explicitly in Eq. (23), the requirement is that the 
change in over a distance of order M(</>) -1 be small.) This local functional takes the 
form of a derivative expansion 



W(<P) = J 



d 4 x 



V(<p) + ±Z(<j ) )(d^) 2 + 



(24) 



Although this is similar in form to the expansion for S'eff, the functions entering the two 
expansions are not the same. In line with the remarks made previously, V((p) and Z{<p) 
differ from Vefj(0, 0) and Z^{cj), 0) by the omission of graphs with internal ^-lines. In 
particular, to 0(g 4 ) we have 



but at the next order 



Vg* ± (Kff) 



(25) 
(26) 



because the second graph in Fig. 1 contributes to V e s but not to V. 

Returning to the bubble nucleation problem, let us define an action 



W (</>) = 



d x 



(27) 



that does display the correct vacuum structure and that, at least for slowly varying 0, is 
a good approximation to W. We can now attempt to evaluate the path integral over </> in 
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Eq. (22) by expanding about the stationary points <j> of Wo. These include the homogeneous 
false vacuum and a bounce solution obeying 



9V 9 4 



(28) 



Defining rj(x) = 4>{x) — <j>(x), we obtain 



W{4>) = W(4>) + d 4 zW'((j);z)r](z) + - d 4 z d 4 z' W" \(f>; z, z')r)(z)r)(z') + 0{rf) (29) 



where primes denote variational derivatives. Note that W'(<f);z) does not vanish, since 
is a stationary point of Wo but not of W. Note also that we cannot use the derivative 
expansion of W inside the path integral, since the integral includes rapidly varying field 
configurations. Inserting Eq. (29) into the path integral and then proceeding as in the 
standard case, we eventually obtain 



e~ Cl e° 2 



-1/2 



J(l + -.-) 



where 



and 



Ci = W{cj> h ) - W(0 fv ) 



c 2 = w'iMiW'iM-'w'M - (<p h ^ fv ). 



(30) 



(31) 



(32) 



Let us now examine the various terms in this expression, beginning with C±. Scaling 
arguments of the type described previously show that the bounce solution to Eq. (28) has 
a characteristic spatial size of order 1 / (g 2 cr) and is therefore slowly varying relative to 
the scale M(<p) entering the one-loop approximation to W, Eq. (23). This both allows us 
to use the derivative expansion (24) to evaluate W((j)b) and also suppresses the higher- 
derivative terms in this expansion. Doing the derivative expansion, one finds that the 
leading contribution to C\ is 



Bn 



1 



d 4 x 



1 



0V> 6 ) 2 + - (0b ^ 0fv) 



1 



(33) 



Recalling that V g 4 = (V e ff) g 4, we see that this is just the result one would expect from 
simply replacing V by V e g in the standard formalism. The next-to-leading contribution is 



B 1 = / d 4 x 



1 „ 

-z g 2((p b )(d^4> b ) 2 + v g 6(4> b ) - (0 b -> 0f v ) 



(34) 



Because V g & ^ {V e ^) g & (although Z g 2 and [Z^) g i are equal in our model), this term is not 
simply the next correction to S e ff. The term beyond this, which is of order unity, involves 
the 0(g 8 ) contributions to V and the 0(g 4 ) contributions to Z, as well as the leading 
four-derivative terms. 
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For C2 we need both W and W" . The derivative expansion can be used for the former; 
because is a stationary point of Wq, the leading contribution here is from 5W = W — Wq. 
Matters are less simple for W"{4>;z, z'). Although the derivative expansion can be used 
when \z — z'\ is large, the behavior for small \z — z'\ is sensitive to the high momentum 
modes and so the derivative expansion fails no matter how slowly varying is. However, 
the relation 

W" = Wq [l + (Wq)~ 1 SW"] (35) 

can be used to obtain formal expansions for (W")~ l and detW" as power series in 
(Wq)~ 1 SW" . The actual utility of these expansions depends on the size of the contri- 
bution from the region of small \z — z'\. In the calculation of C 2 this contribution is 
subdominant and (W") -1 can be approximated by (Wq) -1 ; one finds that C 2 is of order 
unity. For the determinant factor, on the other hand, more terms must be retained: 

det'[W"] = det'[W£] det[/+ (Wq)~ 1 5W"\ 

= det'[Wg] exp {trln[/+ (Wq)~ 1 5W"]} 

= det'TO exp \tY{W^)- l 6W" + ^ tr [(W^)- 1 8W"f + ■■■ 

In the last line, the first term in the exponent is of order 1/g 2 while the second is of order 
unity. In fact, the former term is given by the second graph of Fig. 1 which, it will be 
recalled, contributes to V e s but not to V. It combines with the terms in B\ to give the 
full 0(g 6 ) contribution to the effective potential. Combining all other factors together in 
a prefactor A of order unity, we obtain 

T = Aa A e~( Bo+B ^ (37) 



where 



Bn = Bn 



I 



d 4 x 



- (0b -> 0fv) 



B 1 



(Kff) 3 6(0 b ) + -(Z^ g 2(<f> h )(d^ h )' 



- ((p h -> 0f v ) 



9 Z 



(38) 



(39) 



After seeing how the expansion of the determinant combines with the exponent factors 
to reconstruct S e R to leading and next-to-leading orders, it is natural to speculate that 
this process might continue to higher orders. It does not. In particular, it is at the next 
order that we would encounter the potentially complex contributions to V e s which arise 
from graphs with a single loop. These would appear in a derivative expansion of the 
determinant factor 

-1/2 

(40) 



det'VTOb) 



detW^'( 



>fvj 



However, a derivative expansion of det Wq is valid only for fields which vary slowly relative 
to (V e ff)g 4 . The bounce solution does not satisfy this condition, and so the complex terms 
in V e g do not appear. (One can make a derivative expansion of the determinant in the 
denominator, but (Kff) £ ,s(0f v ) is real.) 
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SCALAR QUANTUM ELECTRODYNAMICS 



These methods can be applied to gauge theories with only minor modifications. To be 
specific, consider the case of scalar electrodynamics with the quartic scalar self-coupling 
taken to be 0(e 4 ) so that the vacuum structure is determined by the one- loop corrections 
to the effective potential? Directly following the approach used for the scalar field example, 
one would integrate out the photon field at the start to obtain an effective action for the 
complex scalar field. The classical equations following from this action have a bounce 
solution of the desired type which gives the leading approximation to V. However, the 
calculation of the next order terms turns out to be much more complicated than in the 
scalar case. These complications can be avoided by an alternative approach. Since the 
bounce solution can be chosen to be entirely real, it is possible to integrate out both 
the photon field and the imaginary part of the scalar field to give an action W{<p) which 
depends only on a single real scalar field. The analysis then proceeds very much as before. 
The 0(e 4 ) terms in V e s and V are identical and, together with the dominant gradient 
term, lead to a contribution to the exponent proportional to e -4 . At order e 6 there are 
two graphs contributing to V e s, neither of which appears in V. Both are recovered from 
the expansion of detW and combine with (Z^) e 2 = (Z) e 2 to give a contribution to the 
exponent which is 0(e~ 2 ). V e s becomes complex at order e 8 but, as before, the failure 
of the derivative expansion of det Wq about the bounce solution prevents the offending 
terms from entering the nucleation rate calculation. In fact, the appearance of infrared 
divergences at small 4> means that the derivative expansion cannot be carried out beyond 
the four-derivative terms. Instead, one must extract the potential and two-derivative terms, 
leaving a remainder which gives an 0(1) contribution to V that can be absorbed in the 
pref actor. 

The issue of gauge-dependence also arises. The nucleation rate is a measurable physi- 
cal quantity and so should be independent of gauge. However, the formulas above express it 
in terms of the effective potential and other gauge-dependent 7-10 quantities. The Nielsen 
identities 8 indicate how this conflict can be resolved. In a class of gauges with photon 
propagator 

JWfc 2 ) = ~%^ +^ (41) 



these identities take the form 

(We do not need the explicit form of the functional C[<f>(x)], but only the fact that it is 
0(e 2 ).) For a uniform field this reduces to 

e^-C(«^. (43, 

Now let us expand this in powers of e 2 to yield a series of identities. The first of these 
states that (V e s) e 4, is gauge-independent, as is easily verified. The next is 

i sjY^ = CM) ajv0 L m 
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Further identities are obtained by making derivative expansions of both sides of Eq. (42) 
and expanding the various terms in powers of e 2 . The first new identity obtained in this 
manner is 



dC e2 



* * (45) 

Let us now apply these identities to V. The leading (0(l/e 4 )) term in the exponent 
involves (V e ff) e 4, and is manifestly gauge-independent. The next term in the exponent, of 
order 1/e 2 , contains (14ff) e e and [Z^) e i. Using the identities we have just obtained, we 
find that the gauge dependence of this term is given by 



d_ 



d 4 X [(Z^) e2 (d^) 2 + (V eS ) e e] 



d 4 x 
d 4 x 



9C e 2 2 d(V eff 



d xC e 2 



dc/> 

(d^C e 2)(d^)+C e 2 

d(V eS ) e * 



d(j) 

d(Ves) t 

dcf> 



(46) 



"□0 + 



dcf> 



Because (pb is a solution of the Euclidean field equations, the last line vanishes. 



CONCLUSION 

In this talk I have described how the decay rate of a metastable vacuum can be 
calculated in a theory whose vacuum structure is determined by radiative corrections. 
As in the standard case, the result may be written as a dimensionful prefactor times 
the exponential of an action involving a bounce solution. To leading approximation this 
exponent is just the tree-level action supplemented by the dominant one-loop contribution 
to the effective potential. The first correction to the exponent arises from the next-to- 
leading contributions to the effective potential and the leading correction to the tree-level 
kinetic part of the effective action. Although smaller than the leading terms, these give 
an addition to the exponent which is larger than order unity and is thus more important 
than the prefactor. It does not appear that this correction need have any particular sign, 
but rather that it might increase the nucleation rate in some theories and reduce the rate 
in others. 

All further corrections may be absorbed into the prefactor. Although some of these 
can be identified with particular terms in the effective action, this is not true of all higher 
corrections. Specifically, the graphs which give rise to complex terms in the effective 
potential cannot, when calculated in the background of the bounce, be expanded in a 
derivative expansion. Consequently, the imaginary part of the effective potential does not 
explicitly enter the bubble nucleation calculation and the problems of interpretation which 
it would entail are avoided. 

This work was supported in part by the U. S. Department of Energy. 
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